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Abstract: 

The complete electroweak 0(a) corrections are calculated for the charged-current 
four-fermion production processes e + e~ — > u T T + fi~u fl , ud/i - ^, and udsc. The calculation 
is performed using complex gauge-boson masses, supplemented by complex couplings to 
restore gauge invariance. The evaluation of the occurring one-loop tensor integrals, which 
include 5- and 6-point functions, requires new techniques. Explicit numerical results are 
presented for total cross sections in the energy range from the W-pair-production threshold 
region up to a scattering energy of 2 TeV. A comparison with the predictions based on 
the "double-pole approximation" (DPA) provided by the generator RacoonWW reveals 
corrections beyond DPA of < 0.5% in the energy range 170—300 GeV, in agreement with 
previous estimates for the intrinsic DPA uncertainty. The difference to the DPA increases 
to 1—2% for y/s ~ 1—2 TeV. At threshold, where the DPA becomes unreliable, the full 
0(a) calculation corrects an improved Born approximation (IBA) by about 1.6%, also 
consistent with an error estimate of the IBA. 
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1 Introduction 



At LEP2, W-pair-mediated four-fermion (4/) production was experimentally explored 
with quite high precision (see Ref. [E] and references therein). The total W-pair cross 
section was measured from threshold up to a centre-of-mass (CM) energy of 207 GeV; 
combining the cross-section measurements a precision of ~ 1% was reached. The W-boson 
mass M\v was determined from the threshold cross section with an error of ~ 200 MeV and 
by reconstructing the W bosons from their decay products within ~ 40 MeV. Deviations 
from the Standard Model (SM) triple gauge-boson couplings, usually quantified in the 
parameters Agf, Ak 7 , and A 7 , were constrained within a few per cent. 

The LEP2 measurements had set the scale in accuracy in the theoretical predictions 
for W-pair-mediated 4/ production. The theoretical treatment and the gained level in 
precision are reviewed in Refs. [ 121 E] • The W bosons are treated as resonances in the full 
4/ processes, e + e~ — > 4/ (+7). Radiative corrections are split into universal and non- 
universal corrections. The former comprise leading-logarithmic corrections from initial- 
state radiation (ISR), higher-order corrections included by using appropriate effective 
couplings, and the Coulomb singularity. These corrections can be combined with the 
lowest-order matrix elements easily. The remaining corrections are called non-universal, 
since they depend on the process under investigation. For LEP2 accuracy, it was suffi- 
cient to include these corrections in the so-called double-pole approximation (DPA), where 
only the leading term in an expansion about the poles in the two W-boson propagators 
is taken into account. Different versions of such a DPA have been used in the literature [ 
Hi El El 13 IE] • Although several Monte Carlo programs exist that include universal correc- 
tions, only two event generators, YFSWW [El El and Racoon WW [013110], include 
non-universal corrections. 

In the DPA approach, the W-pair cross section can be predicted within ~ 0.5% (0.7%) 
in the energy range between 180 GeV (170 GeV) and ~ 500 GeV, which was sufficient for 
the LEP2 accuracy of ~ 1% for energies 170— 207 GeV. In the threshold region (y/s < 
170 GeV), the DPA is not reliable, and the best available prediction results from an 
improved Born approximation (IBA) based on leading universal corrections only, and thus 
possesses an intrinsic uncertainty of ~ 2%. At energies above 500 GeV effects beyond 
0(a), such as Sudakov logarithms at higher orders, become important and should be 
included in predictions at per-cent accuracy. 

At a future International e + e~ Linear Collider (ILC) [ El EH E3] , the accuracy of the 
cross-section measurement will be at the per-mille level, and the precision of the W-mass 
determination is expected to be ~ 10 MeV [E] by direct reconstruction and ~ 7 MeV from 
a threshold scan of the total W-pair-production cross section [EJE]- The theoretical 
uncertainty (TU) for the direct mass reconstruction at LEP2 is estimated to be of the order 
of ~ 5 MeV [H5| to < 10 MeV [EE], based on results of YFSWW and Racoon WW; 
theoretical improvements are, thus, desirable for an ILC. For the cross-section prediction 
at threshold the TU is only ~ 2%, because it is based on an IBA, and thus is definitely 
insufficient for the planned precision measurement of M\y i n a threshold scan. The main 
sensitivity of all observables to anomalous couplings in the triple gauge-boson vertices 
is provided by the W-pair production angle distribution. The TU in constraining the 
parameter A 7 was estimated to be ~ 0.005 [ IT7| for the LEP2 analysis. Since a future ILC 
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is more sensitive to anomalous gauge-boson couplings than LEP2 by more than an order 
of magnitude, a further reduction of the uncertainties resulting from missing radiative 
corrections is necessary. In summary, these considerations demonstrate the necessity of 
a full one-loop calculation for the e + e~ — > 4/ process and of further improvements by 
leading higher-order corrections. 

In this paper we present first results of a complete 0(a) calculation (improved by 
higher-order ISR) for the 4/ final states i/ r r + ju"~P j[i , ud/i~z/^, and udsc, which are relevant 
for W-pair production. 1 The actual calculation is rather complicated. 2 Technically the 
occurring one-loop tensor integrals comprise 5- and 6-point functions up to rank 3, and 
conceptually the W-boson resonances require a treatment in loop diagrams that preserves 
gauge invariance. Details of our approach will be described in a forthcoming publica- 
tion; in the next section we just sketch the most important features of the calculation. 
In Section |H] we present explicit numerical results on total cross sections for scattering 
energies from near the W-pair-production threshold up to 2 TeV. Particular attention is 
paid to a comparison with the DPA and IBA approaches used at LEP2. 

2 Method of calculation 

The actual calculation builds upon the Racoon WW approach [ Ej, where real- 
photonic corrections are based on full matrix elements and virtual corrections are treated 
in DPA. Real and virtual corrections are combined either using two-cutoff phase-space 
slicing or employing the dipole subtraction method [ EOl IS] for photon radiation. We 
also include leading-logarithmic initial-state radiation (ISR) beyond 0(a) in the structure- 
function approach (Ref. and references therein). The presented calculation only differs 
in the treatment of the virtual corrections from RacoonWW. We neglect the masses of 
the fermions whenever possible, i.e. everywhere but in the mass-singular logarithms, and 
set the quark-mixing matrix to the unit matrix. 

In the following we sketch the main difficulties in the calculation and briefly explain 
our solutions. 

2.1 Technical issues 

In contrast to the DPA approach, the one-loop calculation of an e + e~ — ► 4/ pro- 
cess requires the evaluation of 5- and 6-point one-loop tensor integrals. Some 6-point 
diagrams are shown in Figure ^ for illustration. For the generic /1/2/3/4 final state, 
where f\ and ^3 are different fermions excluding electrons and electron neutrinos and /2 
and fi their isospin partners, there are 40 hexagon diagrams, 112 pentagon diagrams, 
and 227 (220) box diagrams in the conventional 't Hooft-Feynman gauge (background- 
field gauge). We calculate the 6-point integrals by directly reducing them to six 5-point 

Electrons and/or positrons in the final state are not yet considered; they deserve further refinements, 
in particular the inclusion of finitc-electron-mass effects in the domain of forward-scattered e ± . We also 
do not yet include final states that can also be produced via two resonant Z bosons, so called mixed 
CC/NC reactions. These can be taken into account in lowest order in RacoonWW. 

2 Some of the problems appearing in a first attempt of such a calculation were already described in 
Ref. [HHj. Recently the authors of the GRACE/1-LOOP system reported on progress towards a full 
one-loop calculation for e + e~ — > /i~/> M ud in Ref. [EH so that one can expect that the system will be able 
to deal with e + e~ — ► 4/ processes at one loop in the near future. 
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Figure 1: Ten 6-point diagrams contributing to e + e~ — > /i 72/3/4- The remaining 30 6- 
point diagrams are obtained by reversing the fermion flow in one or both of the fermion 
chains corresponding to the outgoing fermions. 

functions, as described in Refs. [ E21 E3] ■ The 5-point integrals are reduced to five 4-point 
functions following the method of Ref . [ 121] . These reduction steps involve so-called (mod- 
ified) Cayley determinants, the zeroes of which are related to the Landau singularities of 
the (sub-) diagrams. We did not encounter numerical problems with these determinants. 
Note that this reduction of 5- and 6-point integrals to 4-point integrals does not involve 
inverse Gram determinants composed of external momenta, which naturally occur in the 
Passarino-Veltman reduction of tensor to scalar integrals. The latter procedure leads 
to serious numerical problems when the Gram determinants become small, which happens 
usually near the boundary of phase space but can also occur within phase space because 
of the indefinite Minkowski metric. 

We use, however, Passarino-Veltman reduction to calculate tensor integrals up to 
4-point functions, which involves inverse Gram determinants built from two or three mo- 
menta. This, in fact, leads to numerical instabilities in phase-space regions where these 
Gram determinants become small. For these regions we have worked out two "rescue sys- 
tems" : one makes use of expansions of the tensor coefficients about the limit of vanishing 
Gram determinants; in the other, alternative method we numerically evaluate a specific 
tensor coefficient, the integrand of which is logarithmic in Feynman parametrization, and 
derive the remaining coefficients as well as the scalar integral from it algebraically. This 
reduction again involves only inverse Cayley determinants, but no inverse Gram determi- 
nants. 

In addition to the evaluation of the one-loop integrals, also the evaluation of the 
three spinor chains corresponding to the three external fermion-antifermion pairs is non- 
trivial, because the chains are contracted with each other and/or with four-momenta in 
many different ways. There are O(10 3 ) different chains to calculate, so that an algebraic 
reduction to a standard form which involves only very few standard chains is desirable. 
We have worked out algorithms that reduce all occurring spinor chains to O(10) standard 
structures without introducing coefficients that lead to numerical problems. 

2.2 Conceptual issues 

The description of resonances in (standard) perturbation theory requires a Dyson 
summation of self-energy insertions in the resonant propagator in order to introduce the 
imaginary part provided by the finite decay width into the propagator denominator. It is 
well known that this procedure in general violates gauge invariance, i.e. destroys Slavnov- 
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Taylor or Ward identities and disturbs the cancellation of gauge-parameter dependences, 
because different perturbative orders are mixed (see, for instance, Ref. [0] and refer- 
ences therein). Several solutions have been described for lowest-order predictions, but 
the general problem is still considered as unsolved for a consistent evaluation of radiative 
corrections. The DPA provides a gauge-invariant answer in terms of an expansion about 
the resonance, 3 but in the full calculation we are after a unified description that is valid 
both for resonant and non-resonant regions in phase space, without any matching between 
different treatments for different regions. 

In our calculation we use a generalization of the so-called "complex-mass scheme" 
(CMS), which was introduced in Ref. for lowest-order calculations, to the one- loop 
level. In this approach the W- and Z-boson masses are consistently considered as complex 
quantities, defined as the locations of the poles in the complex p 2 plane of the correspond- 
ing propagators with momentum p. Gauge invariance is preserved if the complex masses 
are introduced everywhere in the Feynman rules, in particular, in the definition of the 
weak mixing angle, 



which is now derived from the ratio of the complex mass squares. The (algebraic) rela- 
tions, such as Ward identities, that follow from gauge invariance remain valid, because 
the gauge-boson masses are modified only by an analytic continuation. Since this contin- 
uation already modifies the lowest-order predictions by changing the gauge-boson masses, 
double-counting of higher-order effects has to be carefully avoided by an appropriate 
renormalization procedure. 

The use of complex gauge-boson masses necessitates the consistent use of these com- 
plex masses also in loop integrals. To this end, we have derived all relevant one-loop inte- 
grals with complex internal masses. The IR-singular integrals were taken from Ref. [177]. 
Concerning the non-IR singular cases, we have analytically continued the results of 
Ref. [ 28j for the 2-point and 3-point functions 4 , and the relevant results of Ref. [|2H] for 
the 4-point functions. We have checked all these results by independent direct calculation 
of the Feynman-parameter integrals. 

2.3 Checks 

The complexity of the calculation enforces a number of consistency checks to prove 
the reliability of the results. We have performed the following checks: 

• UV finiteness is checked by verifying the independence of the (arbitrary) reference 
scale \x of dimensional regularization. 

• IR finiteness is checked by varying the logarithm In A of the (formally infinitesimal) 
photon mass A, which leaves the sum of the virtual and the soft-photonic corrections 

3 The recently proposed approach to describe unstable particles within an effective field theory is 
equivalent to a pole expansion. 

4 Note that the result of Ref. [EHj for the scalar two-point function is not valid in general for complex 
masses. In this case an extra rj function has to be added. The same comment applies to the results for 
the 2-point tensor integrals in Ref. [125). 
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(slicing approach) or of the virtual and endpoint contributions (i.e. the singular 
parts that are subtracted from the virtual corrections in the subtraction approach) 
invariant. 

• Mass singularities related to collinear photon emission or exchange are checked by 
verifying the independence of the sum of the virtual corrections and the subtraction 
endpoint contributions from the small masses of the external fermions. 

• Gauge invariance is checked by comparing the result obtained within the 't Hooft- 
Feynman gauge with an independent result obtained within the background-field 
formalism [|3U]. Apart from diagrams that involve only fermion-gauge-boson cou- 
plings in the loops, the contributions of individual Feynman graphs are in general 
different in the two approaches. 

• The real corrections are taken over from RacoonWW [0 EJ> where they were 
checked by two independent calculations in detail. Moreover, agreement was found 
between the results obtained with phase-space slicing or dipole subtraction. 

• The scalar loop integrals for complex masses have been calculated in two completely 
independent ways and implemented in two independent in-house libraries. We 
checked that these results agree with each other and in the limit of zero width 
also with FF [EH- 

• Two completely independent calculations have been performed within our group, 
revealing good agreement. All algebraic manipulations, including the generation of 
Feynman diagrams and the reduction of amplitudes to standard forms, have been 
done using two independent programs. 5 The evaluations of all scalar and tensor loop 
integrals are based on two independent in-house libraries, which employ the two 
different rescue systems mentioned above. We consider this as the most important 
and convincing check. 

We emphasize that all these checks, including the gauge-invariance check, have been 
carried out for non-zero gauge-boson widths. 

5 The amplitudes are generated with FeynArts, using the two independent versions 1 and 3, as 
described in Refs. [ (321 an d [ESI; respectively. The algebraic manipulations are performed using two in- 
dependent in-house programs implemented in Mathematica, one of which builds upon FormCalc [ IH4| . 
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3 Numerical results 



3.1 Input parameters and setup 

For the numerical evaluation we use the following set of SM parameters, 





= 1.16637 x 10- 5 GeV" 2 , 


a(0) 


= 1/137.03599911, 




= 0.1187, 


M w 


= 80.425 GeV, 


M z 


= 91.1876 GeV, 


r z 


= 2.4952 GeV, 


M H 


= 115 GeV, 










m e 


= 0.51099892 MeV, 


nip 


= 105.658369 MeV, 


m T 


= 1.77699 GeV, 


m u 


= 66 MeV, 


m c 


= 1.2 GeV, 


m t 


= 178 GeV, 


m d 


= 66 MeV, 


m s 


= 150 MeV, 


m h 


= 4.3 GeV, 



(3.1) 

which essentially follows Ref. [[32]- For the top-quark mass m t we have taken the more 
recent value of Ref. [EE]- The masses of the light quarks are adjusted to reproduce 
the hadronic contribution to the photonic vacuum polarization of Ref. [EZ1- Since we 
parametrize the lowest-order cross section with the Fermi constant G M (G M scheme), i.e. 
we derive the electromagnetic coupling a according to = \flG ^M^iX — M^/MD/tt, 
the results are practically independent of the masses of the light quarks. Moreover, 
this procedure absorbs the corrections proportional to mf/M^j in the fermion-W-boson 
couplings and the running of ot(Q 2 ) from Q 2 = to the electroweak scale. In the relative 
radiative corrections, we use, however, a(0) as coupling parameter, which is the correct 
effective coupling for real photon emission. 

QCD corrections are treated in the "naive" approach of multiplying cross sections 
and partial decay rates by factors (1 + a s /vr) per hadronically decaying W boson. The 
W-boson width T-w is calculated from the above input including electroweak 0(a) and 
QCD corrections, yielding 

T w = 2.09269848 . . . GeV. (3.2) 

For the full one-loop calculation and for the DPA approach this procedure ensures that 
the effective branching ratios for the leptonic, semileptonic, and hadronic W decays, which 
result from the integration over the decay fermions, add up to 1. In order to keep this 
normalization also for the IBA, T-w is calculated in the corresponding approximation, i.e. 
in the G M scheme without electroweak corrections, yielding 



w 



IBA 



2.10009936... GeV. (3.3) 



A detailed description of the IBA, as used in Racoon WW, can be found in Ref. [E 

In the following we discuss only total cross sections without any phase-space cuts. The 
presented results have been obtained with 10 7 events, using the subtraction method. 



3.2 Results for total cross sections 

Tables show some representative results on total cross sections for the final states 
v t t^ n~v^ ud/i - ^, and udsc in various approximations for different CM energies a/s. The 
numbers in parentheses represent the uncertainties in the last digits of the predictions. 
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e + e — > u T T + jj, v i 



v^/GeV 


Born(FW) 


Born(CMS) 


IBA 


DPA 


ee4f 


161 


50 04(2) 


50 01(2) 

\J KJ • KJ -L \ J 


37 18(2) 


37 08(2) 

-.7 1 . V7 v7 \ / 


37 95(2) 

"7 1 . t_7 ^7 \ J 






[-0.06%] 

L J 


[-25.67(6)%] 


[-25.90(3)%] 

L \ / J 


[-24.12(4)%] 


170 


160 53(6) 


160 44(6) 


129 12(6) 


129 17(6) 


129 23(6) 






[-0.06%] 


[-19.52(5)%] 


[-19.53(3)%] 

L \ / J 


[-19.45(3)%] 

L \ / J 


189 


216 57(8) 

-L V7 • .7 1 \ V_7 I 


216 45(8) 


191 89(8) 

-L t_7 -L . ^7 t_7 \ ^7 / 


191 66(9) 


190 89(9) 

-L ' 7 V 7 • v7 '.7 \ '.7 J 






[-0 06% 


[-11 35(5)%1 


[-11 50(2)%1 


[-11 81(3)%] 


200 


220.41(9) 


220.29(9) 


201.13(9) 


200.04(10) 


199.21(10) 






[-0.06%] 


[-8.70(6)%] 


[-9.24(2)%] 


[-9.57(3)%] 


500 


86.95(5) 


86.90(5) 


92.79(5) 


89.81(6) 


89.13(6) 






[-0.06%] 


[+6.78(9)%] 


[+3.29(3)%] 


[+2.57(4)%] 


1000 


33.35(2) 


33.33(2) 


38.04(4) 


35.76(3) 


35.37(3) 






[-0.06%] 


[+14.12(14)%] 


[+7.21(5)%] 


[+6.12(6)%] 



Table 1: Total cross sections in fb for e + e~ — > v t t + in Born approximation (in the 
fixed-width and complex-mass schemes), IBA, DPA, and using the full 0(a) correction 
(ee4f); all but the Born cross sections include higher-order ISR corrections. 



e + e — > ud/i u l 



v/i/GeV 


Born(FW) 


Born(CMS) 


IBA 


DPA 


ee4f 


161 


150.15(7) 


150.07(7) 


115.75(7) 


115.48(7) 


118.12(8) 






[-0.06%] 


[-22.87(6)%] 


[-23.09(3)%] 


[-21.29(4)%] 


170 


481.6(2) 


481.3(2) 


402.0(2) 


401.8(2) 


402.1(2) 






[-0.06%] 


[-16.48(5)%] 


[-16.58(3)%] 


[-16.47(3)%] 


189 


649.7(3) 


649.4(3) 


597.4(3) 


596.1(3) 


593.8(3) 






[-0.06%] 


[-8.00(5)%] 


[-8.26(3)%] 


[-8.55(3)%] 


200 


661.3(3) 


660.9(3) 


626.2(3) 


622.2(3) 


619.7(3) 






[-0.06%] 


[-5.26(6)%] 


[-5.91(3)%] 


[-6.24(3)%] 


500 


260.9(1) 


260.8(1) 


288.9(2) 


279.6(2) 


277.4(2) 






[-0.06%] 


[+10.78(9)%] 


[+7.14(4)%] 


[+6.37(4)%] 


1000 


100.10(6) 


100.04(6) 


118.44(13) 


111.45(9) 


110.17(10) 






[-0.06%] 


[+18.39(15)%] 


[+11.34(5)%] 


[+10.13(7)%] 



Table 2: Total cross sections in fb for e + e~ — > ud/i~z/ M in Born approximation (in the 
fixed-width and complex-mass schemes), IBA, DPA, and using the full 0(a) correction 
(ee4f); all but the Born cross sections include higher-order ISR and QCD corrections. 
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e + e — > udsc 



v^/GeV 

V / 


Born(FW) 


Born(CMS) 


IBA 


DPA 


ee4f 


161 


450.5(2) 


450.3(2) 


360.4(2) 


359.7(2) 


367.6(3) 






[-0.06%] 


[-19.97(6)%] 

L \ / J 


[-20.16(4)%] 


[-18.37(4)%] 


170 


1444.9(5) 


1444.1(5) 


1251.6(5) 


1250.1(6) 


1250.7(6) 






[-0.06%] 

L J 


[-13.33(5)%] 


[-13.49(3)%] 


[-13.39(3)%] 


189 


1949.3(8) 


1948.2(8) 


1860.0(8) 


1853.8(9) 


1846.9(9) 






[-0.06%] 

L J 


[-4.53(6)%] 

L \ / J 


[-4.90(3)%] 


[-5.20(3)%] 

L \ / J 


900 


1 08^ 

IZ/OO. Z/\Oj 


1 080 Q(Q\ 
-L c/O^.iy I O J 


1 Q4Q ^(Q) 




1 Q97 7(Q} 






[-0.06%] 


[-1.68(6)%] 


[-2.45(3)%] 


[-2.78(3)%] 


500 


782.9(4) 


782.5(4) 


899.4(5) 


869.4(6) 


862.9(6) 






[-0.06%] 


[+14.94(9)%] 


[+11.05(4)%] 


[+10.27(4)%] 


1000 


300.3(2) 


300.1(2) 


368.7(4) 


346.1(3) 


342.1(3) 






[-0.06%] 


[+22.86(16)%] 


[+15.26(5)%] 


[+13.98(7)%] 



Table 3: As in Tabled but for the process e + e — > udsc. 

Columns two and three in each table contain the two versions of the lowest-order cross 
section for the full e + e~ — > 4/ processes corresponding to the different treatments of finite- 
width effects as provided by the "fixed- width scheme" (FW) and the complex-mass scheme 
(CMS). In the FW scheme the finite constant width, and thus the complex mass, is only 
inserted into the resonant propagators. The relative difference 0"B O rn(CMS) / crBom(FW) — 1 
of the schemes in lowest order is given by the numbers in square brackets in the third 
columns. Note that the difference is only 0.06% for the considered energies, so that it is 
not essential to which Born cross section we normalize relative corrections. We have not 
given an error on this difference, because the two Born predictions are strongly correlated. 
The fourth columns show the IBA [ [HI] implemented in RacoonWW, which is based 
on universal corrections only and includes solely the contributions of the CC03 diagrams; 
the numbers in square brackets are defined as <5iba = criBA/o"Born(CMS) — 1. The fifth 
columns show the DPA of RacoonWW 6 , which comprises also non-universal corrections 
[E|; the numbers in square brackets are defined as <5dpa = cr DPA/c r Born(FW) — 1. We 
normalize odpa to (T Bom (FW), because the lowest-order part of the DPA is per default 
evaluated in the FW scheme in RacoonWW. Finally, the last columns (ee4f) contain 
the full one- loop corrections to e + e~ — > 4/; the numbers in square brackets are defined 
as 5 ee 4f = cr ee 4f/crBora(CMS) — 1. Here we normalize to <7Born(CMS), because the full 
e + e~ — > 4/ calculation is consistently performed in the CMS. Note that all but the "Born" 
numbers include improvements by ISR effects beyond 0(a), as described in Ref. [Ej. 
Additionally the results on the semileptonic and hadronic cross sections (all but the Born 
cross results) include naive QCD corrections, as explained in Section I3~T1 For better 

6 We recall that the DPA of RacoonWW goes beyond a pure pole approximation in two respects. 
The real-photonic corrections are based on the full e + e _ — > 4/ + 7 matrix elements, and the Coulomb 
singularity is included for off-shell W bosons. Further details can be found in Ref. [Ej. 
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illustration Figure El depicts the predictions for the energy ranges of LEP2 and of the 
high-energy phase of a future ILC, focusing on the leptonic final state v t t + The 
respective figures for the relative corrections 5 to the semileptonic and hadronic final 
states look almost identical, up to an offset resulting from the QCD corrections. 

A comparison between the DPA and the predictions based on the full 0(a) corrections 
reveals differences in the relative corrections 5 of < 0.5% (0.7%) for CM energies ranging 
from y/s ~ 170 GeV to 300 GeV (500 GeV). This is in agreement with the expected 
reliability of the DPA, as discussed in Refs. [EHEIIZI- At higher energies, the deviations 
increase and reach 1—2% at yfs = 1— 2TeV. In the threshold region (y/s < 170 GeV), as 
expected, the DPA also becomes worse w.r.t. the full one-loop calculation, because the 
naive error estimate of (cx/ir) x (r^/Mw) times some numerical safety factor of (9(1 — 10) 
for the corrections missing in the DPA has to be replaced by (a/n) x T w / (-y/s — 2Mw) in 
the threshold region and thus becomes large. In view of that, the DPA is even surprisingly 
good near threshold. For CM energies below 170 GeV the LEP2 cross section analysis was 
based on approximations like the shown IBA, which follows the full one-loop corrections 
even below the threshold at y/s = 2M\y within an accuracy of about 2%, as expected in 
Ref. [Unj. The shape of the relative corrections in the threshold region is determined by 
ISR. The minimum in the relative corrections is correlated to the maximum in the slope 
of the total cross section. 

A detailed comparison [ |3| between the DPA versions of Ref. [@j, of YFSWW [ 
El , and of RacoonW W [ HH1 showed differences of the same order (or smaller) 
as the differences between the DPA of RacoonWW and the result based on the full 
one-loop calculation of e + e~ — > 4/ discussed above. Therefore, these DPA predictions 
are consistent with the results based on the full one-loop calculation. The DPA version 
presented in Ref. however, deviates from the full one-loop calculation by about 0.9% 
for typical LEP2 energies. 

In addition to the above results which include ISR beyond 0(a) and QCD corrections, 
we present also explicit numbers that include only the genuine 0(a) corrections to facili- 
tate a comparison to calculations made by other groups in the future. Table 0] shows the 
relative 0(a) corrections, i.e. without higher-order ISR improvements and without QCD 
corrections, both for the DPA of RacoonWW and for the full e + e~ — > 4/ calculation 
for the three final states zvr + /i~z/ M , ud/i~z/^, and udsc. 

3.3 Remaining theoretical uncertainties 

We have reduced the TU for the charged- current processes e + e~ — > v T r + pTv^ xsApTv^ 
udsc, in particular in the threshold region of W-pair production, considerably by calcu- 
lating the full 0(a) corrections. ISR beyond 0(a) is included via structure functions in 
leading-logarithmic accuracy. For energies below ~ 500 GeV, the remaining uncertainties 
resulting from missing electroweak corrections are then dominated by the next-to-leading 
logarithmic electromagnetic corrections of order (a/ir) 2 log(m%/s) which can be estimated 
to contribute < 0.1%. Near the W-pair-production threshold, higher-order effects of the 
Coulomb singularity are still missing. These are estimated to ~ 0.2% [EH1EH]- Thus, we 
estimate the theoretical uncertainty due to unknown electroweak higher-order effects in 
the present calculation to be a few 0.1% from the threshold region to about ~ 500 GeV. 
At higher energies leading and subleading electroweak high-energy logarithms, such as 
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Figure 2: Absolute cross section a (upper plots) and relative corrections 8 (lower plots), as 
defined in the text, to the total cross section without cuts for e + e~ — > u t t + ji^u^ obtained 
from the IBA, DPA, and the full 0(a) calculation (ee4f). All predictions are improved 
by higher-order ISR. 
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1 

e^e — > 


1 




ud/x" 




udsc 


v^/GeV 


W[%] 


5ee4f[%] 


5dpa[%] 


5ee4f[%] 


5dpa[%] 




161 


-31.30(3) 


-29.53(3) 


-31.30(3) 


-29.56(3) 


-31.27(3) 


-29.62(3) 


170 


-21.67(2) 


-21.59(2) 


-21.75(2) 


-21.65(2) 


-21.81(2) 


-21.72(2) 


189 


-11.74(2) 


-12.05(2) 


-11.84(2) 


-12.13(2) 


-11.94(2) 


-12.22(2) 


200 


-9.15(2) 


-9.48(2) 


-9.25(2) 


-9.57(2) 


-9.33(2) 


-9.64(2) 


500 


+3.55(3) 


+2.84(4) 


+3.50(3) 


+2.76(4) 


+3.37(3) 


+2.63(4) 


1000 


+7.08(5) 


+5.98(6) 


+7.16(5) 


+5.98(7) 


+6.89(5) 


+5.69(7) 



Table 4: Genuine relative 0(a) corrections, i.e. without higher-order ISR improvements 
and QCD corrections. 

Sudakov logarithms, beyond one loop have to be taken into account in addition to match 
this accuracy. 

For a thorough estimate of the total theoretical uncertainty an inclusion of QCD 
effects is indispensable for the processes involving final-state quarks. In order to reach 
a precision of the order of a few 0.1% there, it is certainly necessary to improve the 
treatment of 0(a s ) corrections (and beyond), including a proper matching with parton 
showers. Bose-Einstein and colour interconnection effects may also play an important role. 

4 Summary 

We have presented results on total cross sections for the charged-current four-fermion 
production processes e + e~ — > v T T + fj l ~v fl , ud/z~z/ M , udsc which, for the first time, include 
the complete electroweak 0(a) corrections. The calculation is consistently performed 
using complex gauge-boson masses, supplemented by complex couplings to restore gauge 
invariance. The evaluation of the occurring one-loop tensor integrals, which include 5- and 
6-point functions, required new techniques. Moreover, we have developed algorithms to 
reduce the large number of different spinor chains to a set of very few standard structures. 

A comparison with the predictions for the total cross section without cuts based on 
the DPA provided by the generator RacoonWW reveals corrections beyond DPA of 
< 0.5% (0.7%) for CM energies ranging from ^ ~ 170 GeV to 300 GeV (500 GeV), 
consistent with previous estimates on the intrinsic DPA uncertainty. The difference to 
the DPA increases to 1—2% for y/s ~ 1— 2TeV. At threshold, where predictions had 
to be based on an IBA at LEP2, the full 0(ct) calculation corrects the IBA by about 
1.6%, also consistent with a previous error estimate. The full 0(a) calculation, improved 
by higher-order effects from ISR, reduces the remaining TU due to unknown electroweak 
higher-order effects to a few 0.1% for scattering energies from the threshold region up to 
~ 500 GeV; above this energy leading high-energy logarithms, such as Sudakov logarithms, 
beyond one loop have to be taken into account to match this accuracy. At this level of 
accuracy, also improvements in the treatment of QCD corrections to semileptonic and 
hadronic e + e~ — > 4/ processes will be necessary in the future. 
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Results for differential cross sections as well as details of the calculation will be pre- 
sented elsewhere. 
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